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ABSTRACT 


Nowadays, inverters are power electronics devices used in the 
renewable energy conversion systems. Therefore, the continuity of 
service of power converter has been a major concern. In this paper, a 
technique for diagnosing multiple open switch faults in three phase 


Accepted Jul 10, 2020 voltage inverters feeding Doubly Fed Induction Generator DFIG is 


presented. It is based on the so-called the Lissajous curves of DFIG 


Keywords: rotor currents. For this purpose, the modeling of a non-defected wind 
DFIG conversion system based on mathematic model created in Matlab 

Simulink is firstly presented, after that, an indirect stator field vector- 
Diagnosis 


oriented control is applied to obtain the wind system performance. 
Finally, rotor currents Lissajous curves are analyzed in case of a non- 
defected and defected inverter. Simulation results show the 
effectiveness of the proposed method. 
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1. INTRODUCTION 

Faced with a growing demand for energy, the industrial world is moving more and more towards 
renewable energies and particularly wind turbines to produce electricity. Currently, wind turbines based on 
the Doubly Fed Induction Generator (DFIG) are the most widely used in terrestrial wind farms. The wind 
energy conversion chain is a device that transforms part of the wind kinetic energy into mechanical energy 
available on a transmission shaft and then into electrical energy via a generator [1]. These installations are 
nowadays in full expansion and are of very considerable importance, whether from the point of view of 
electrical or economical energy requirements. However, these wind turbines may be subject to partial or total 
damage, an untimely shutdown of a wind chain inevitably causes inestimable losses. So, to prevent such a 
shutdown from being caused by electrical or mechanical failures, it is important to detect these defects early. 

There have been many studies on wind energy conversion system. Whether at the level of 
regulation [2], of automatic generation control [3] or of faults detection [4]. Our work is part of a project that 
consists in modeling and simulation of the variable speed wind energy conversion chain based on the DFIG 
in case of a non-defected and defected converter. A lot of works has been done to detect static converters 
faults where many methods were used, of which we can specify: the Neural Network Approach [5], the 
Fuzzy Logic Approach [6] and Multidimensional Data Analysis [7]. In general, those methods are based on 
the trajectory of the current vector (instantaneous or average) using analytical models or pattern recognition 
methods to detect faults in the control circuits or power semiconductor faults (Open/ Short circuit faults 
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affecting the IGBT switches) [7]. The Lissajous curves of rotor currents can also detect various inverter faults 
and this without disturbing the generator operation. These curves were presented in different studies such as 
the rotating machinery field [8], and faults diagnosis [9-10], it was also studied in [11] to detect wind turbine 
doubly fed induction machine under grid defects. The aim of this paper is to present Lissajous curves as a 
diagnosis method for detecting three phase voltage inverters feeding DFIG in the case of multiple open 
switches faults. 


2. CONCEPT OF WIND CONVERSION SYSTEM 

A wind system based on a DFIG is generally given under the topology illustrated by the Figure 1. 
The wind turbine drives the DFIG via a gearbox, where it is connected directly to the grid by the stator and 
by the rotor through two three-phase statics converters controlled by pulse width modulation (PWM). One in 
the DFIG rotor side called Rotor Side Converter (RSC) and the other in the grid side called Grid Side 
Converter (GSC). Thanks to the bidirectionality of these two converters, the DFIG can operate in sub 
synchronous and super synchronous mode. PI type power regulators are used. They allow reaching at the 
same time a good dynamic and a good robustness. Rotor side converter acts as a PWM rectifier during the 
machine working in super synchronous mode and as an inverter during sub synchronous mode [13]. The 
purpose of the RSC is providing decoupled active and reactive stator power regulation Ps and Qs depending 
on the reference torque delivered by the PI controllers. Grid side converter acts as a PWM inverter during the 
machine working in super synchronous mode and as a rectifier during sub synchronous mode [13]. The 
purpose of the GSC is controlling the exchange of power flows with the grid via the rotor, maintaining the 
DC bus at a constant voltage level and imposing the reactive power Qs at 0. 





Figure 1. Wind energy conversion system using a doubly fed induction generator [12] 


3. WIND CONVERSION SYSTEM MODELING 
3.1. Wind turbine modeling 

Our model is based on the steady state power characteristics of the turbine. The rigidity of the drive 
shaft is assumed infinite, the friction coefficient and the inertia of the turbine must be combined with those of 
the generator coupled to the turbine [14]. 
— The Mechanical power recovered by the turbine 





Prurbine= * CÀ) + p ET R3 * p? (1) 
and à = Ee (2) 
where 


Qı Speed of rotation before multiplier 
R : Radius of the aerogenerator 
p : Density of air, 1.225kg/m? 


The coefficient of power Cp represents the wind turbine aerodynamic efficiency and depends on the 
turbine characteristic. Figure 2 represents the variation of this coefficient as a function of the speed ratio and 
the angle of the orientation of the blade B. 

— The aerodynamic torque is: 


P . 
Caero ~~ ara (3) 
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— The fundamental equation of dynamics 


dQ J dwy 
aF Se = Cr — Com — Cris (4) 


= C= Cay Sf Ab 


Where 
J : Total inertia 
0.5 : Speed of rotation after multiplier 
Comi : Electromagnetic torque 
Ce. : Viscous friction torque 
f : Coefficient of viscous friction* 


*The resistant torque due to friction is modeled by a coefficient of viscous friction f [16]. 


The fundamental equation of dynamics allows determining the evolution of the mechanical speed 
from the total mechanical torque applied to the rotor. 





Figure 2. Characteristic of the power coefficient based on lambda variation [15] 


3.2. Modeling of doubly-fed induction machine 
Model of the asynchronous machine in Park reference frame (d, q) are [17]: 
d 
Vsa = Rs * lsa + a Psa) — Ws * Psq 
d 
Vsq = R; * Isq + = (Psa) — Ws * Psa 


d 
Vra = Rr * lya + dt (Pra) — (ws-wr) * Prq 


d (5) 
Vig = R, * Lrq + Z (Pra) — (ws-wr) * Ord 
— The electromagnetic torque equation is given by: 
Cee 
Cem = S p * A (Pra = leg — Ọrq * sa) (6) 


3.3. Field oriented control 
We choose a two-phase reference frame (d, q) linked to the rotating field. The stator flux is oriented 


along the direct axis. So we can write: 


Pas = Ps 
Pas = 0 
M ; 
Cem = ZP. 7, Pas: lgr (7) 
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Active and reactive powers can be written as follows: 








— Vs*M : 
P, =g Ls $ lar 
Vs*M . (8) 
S 


4. LISSAJOUS CURVES FOR DFIG DEFECTS DIAGNOSIS 

Lissajous curves or figures, also called Bowdith curves, are the trajectories of a point on the plane 
whose coordinates each have a sinusoidal movement of period, amplitude and phase which may be different 
[18]. So by expressing a sine wave as a function of another sine wave, interesting patterns emerge. This 
family of curves was studied first by the American mathematician Nathaniel Bowdith in 1815, then in more 
detail by the French mathematician Jules-Antoine Lissajous in 1857-1858. Lissajous used a narrow stream of 
sand pouring from the base of a compound pendulum to produce the curves [19-20]. 

Lissajous curves are described by the parametric equations below [21]: 


x = Asin(at + a) y = Bsin(bt) (9) 


If the two waves are exactly the same frequency (a=b), have the same amplitude (A=B), but out of 
synch by 2/2, then one of the waves will be like a cosine, and the resultant output will be a circle: 


x = Acos(at) y = Asin(at) (10) 
ety a At (11) 


These curves are applied in this paper to the DFIG rotor three-phase currents. For this purpose, it is 
necessary to apply the Concordia transformation [22] to the rotor basic currents Iapc, 1n order to obtain the 
rotor currents Concordia components (I-a, Ig ) such as: 


haj a|! a Ta He 12 
Lp T3 0 v3 _ v3 o ( ) 
2 2 rc 


5. SIMULATION RESULTS 
In order to show the effectiveness of the proposed method, simulations of the wind energy 
conversion chain were carried out on Matlab/Simulink. 


5.1. Free-deffects operation 

This case is a reference to which the defected inverter feeding the DFIG cases are compared. 
Figure 3 represent respectively the system characteristics and rotor currents Lisssajous curve during 
free-defects operation. 


5.2. Discussion 

During the transient regime, the rotor speed is slightly pulsating, presenting beats at the first 
moment of starting, after that the speed reaches the steady state by stabilizing at the nominal value. The curve 
of the electromechanical torque presents at startup large beats in a short time interval then stabilizes at a 
negative value which reflects the generator mode. At startup, the bus voltage increases exponentially for a 
few seconds and then returns to a fixed reference value. The reactive power set point will be maintained at 
zero to ensure a unit power factor on the stator side so as to optimize the quality of the energy returned on the 
grid [23]. The active power set point must keep the power factor of the wind turbine optimal. The Lissajous 
curve of Doubly Fed Induction Generator rotor currents Concordia components, as shown in Figure 3 (b), has 
a shape of two circles of which the one inside has six knots: the number of the inverter switches. 
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Figure 3. Free-defects operation, (a) The dynamic behavior of the DFIG, (b) Lissajous curve 
of DFIG rotor currents 


5.3. Faulty operation 

In this application, we are interested in the hypo-synchronous generator operation, so as to treat only 
the transfer of rotor power in a single direction (Grid to Rotor). The grid side converter will therefore be a 
simple unidirectional rectifier. The DC bus voltage is therefore kept fixed and we are interested only in the 
control of the rotor side converter. A converter such as the PWM inverter is likely to have structural defects 
such as power semiconductor faults (switches) short circuit and/or open circuit [24], and therefore this type 
of malfunction may cause damage to the entire production system if the personnel is not notified and an 
untimely shutdown is produced [25]. A study was devoted to the detection and diagnosis of faults in static 
converters in an aspect closer to the physics of the electronic power component and the module which 
described the mechanisms of the main failures. This study has shown that the temperature of the environment 
as well as temperature variations due to thermal cycles are the main causes of failures [7]. In this paper, 
inverter faulty operations are considered in the case of multiple open switches that occur at the same leg and 
at different legs. 









GRID 
f 
POWER j IVER 
p, L 
Ae pT 
if 
GEAR BOX DFIG | 





TON 


Figure 4. The six switches {S} to S,} of the rotor side converter concerned by our faulty operation study [26] 


Open-circuit faults attacking the IGBT switches may occur when the IGBT is disconnected, is 
damaged, or had a problem in its grid control signal. This type of fault is very difficult to recognize directly 
because the motor can continue to run but with a degradation of its performance due to the occurrence of 
fluctuations in the mechanical parameters as well as a variation in the currents where the currents of the other 
two healthy legs take high values to maintain the average torque and the speed. In the presence of this type of 
fault, the starting of the motor cannot always be guaranteed, everything will depend on the value of the 
torque which can be close to zero for certain rotor position [27]. An IGBT open-circuit fault in a switch is 
introduced by removing its gate signal [28]. The rotor currents Concordia components Lissajous curves 
indicate the inverter faulty situation. A free-defects operation is represented by a six knots circle inside a 
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perfect one. If an inverter fault takes place, there is a change in the rotor currents. Therefore, a change in the 
circle shape proves the presence of inverter defect. 

Inverter Single Switch Fault. The figures below show the rotor currents Lissajous curve in the case 
of different level single switch fault. When the fault occurs, the circle with six knots becomes a circle with 5 
knots Figure 4. The orientation of the circle missed knots depends on the defected switch location. In the case 
of one upper switch fault, the missed knot appears at the bottom (a). In the case of one lower switch fault, the 
missed knot appears at the top (b). 
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Figure 4. Lissajous curves of rotor currents Concordia Components, 
(a) Case of a DFIG fed by S4 switch defected inverter, (b) Case of a DFIG fed by S, switch defected inverter 


Inverter Double Switches Fault in the Same Leg. The Figure 5 shows the rotor currents Lissajous 
curve in the case of the same leg double switches faults. When the fault occurs, the circle with six knots 
becomes a circle with 4 knots, it’s close to a shape of butterfly with knot inside Figure 5. The orientation of 
the circle depends on the defected switches location. 

Inverter Double Switches Fault in different Legs. The Figure 6 shows the rotor currents Lissajous 
curve in the case of different leg double switches fault. When the faults occur, the circle with six knots 
becomes a circle with 4 knots. The orientation of the circle missed knots depends on the defected switches 
location. In the case of two upper switches fault, the missed knots appear at the top (a). In the case of two 
lower switches fault, the missed knots appear at the bottom (b). In the case of one upper switch and one lower 
switch faults, the obtained shape is presented in (c). 
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Figure 5. Lissajous curves of rotor currents, 
(a) Case of a DFIG fed by S4 and S, switches defected inverter, (b) Case of a DFIG fed by S3 and S, 
switches defected inverter, (c) Case of a DFIG fed by S; and Sẹ switches defected inverter 
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Figure 6. Lissajous curve of rotor currents, (a) Case of a DFIG fed byS, and S} switches defected inverter, 
(b) Case of a DFIG fed by S, and S, switches defected inverter, (c) Case of a DFIG fed by S4 and S,switches 
defected inverter 


6. CONCLUSION 

In the first part of this paper we studied the operating mode of the wind energy conversion system 
namely the model of the variable speed turbine and the model of the Doubly Fed Induction Generator. In the 
second part, a method was developed for detecting three phase inverter Single and Multiple Open Switches 
faults, based on the rotor currents Concordia components Lissajous curves. In the case of a free-defects 
operation, the Lissajous curve has a shape of two circles of which the one inside has six knots (exactly the 
number of the inverter switches). The knots number decreases and the circle orientation changes according to 
the number of the inverter defected switches (single or double) and their location. We are still pursuing 
simulations of the inverter feeding Doubly Fed Induction Generator in the case of Short Circuit Switches 
Fault before moving to the Grid Side Converter to do the complete inverter faulty operation study. To 
complete the project, a system of supervision and decision making regarding the type of defect and its 
location must be established. 
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